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Wo consider turbulence in the Gross-Pitacvsky model and study the creation of a coherent con- 
densate via an inverse cascade originated at small scales. The growth of the condensate leads to a 
spontaneous breakdown of symmetries of small-scale over-condensate fluctuations: first, statistical 
isotropy is broken, then series of phase transitions mark the change of symmetry from the two- fold to 
three-fold to four-fold. At the highest condensate level reached, we observe a short-range positional 
and long-range orientational order (similar to a hexatic phase in the condensed matter physics). In 
other words, the longer one pumps the system the more ordered it becomes. We show that these 
phase transitions happen when the driving term corresponds to an instability (i.e. it is multiplica- 
tive in the k-space) but not when the system is pumped by a random force. Thus we demonstrate 
for the first time non-universality of the inverse-cascade turbulence. We also describe anisotropic 
spectral flux flows in k-space, anomalous correlations of fluctuations and collective oscillations of 
turbulence-condensate system. 

PACS numbers: 47.27.Gs; 03.75.Hh; 42.65.Sf 



I. INTRODUCTION 

Probably the most unexpected discovery made in the 
study of turbulence is an inverse cascade. Flying at the 
face of an intuitive picture of turbulence as the process 
of fragmentation and dissipation, inverse cascade moans 
self-organization, i.e. appearance of large-scale motions 
out of a small-scale noise. The process of inverse cas- 
cade culminates in the creation of spectral condensate 
i.e. a mode coherent across the whole system. At this 
point turbulence shares many properties with quantum 
physics as displaying both fluctuations and coherence. 
This closeness shows perhaps most vividly in the studies 
on non-equilibrium states within the framework of Non- 
linear Schrodinger (or Gross-Pitaevsky) Equation, called 
here NSE. In the realm of classical physics, this model de- 
scribes spectrally narrow distribution of nonlinear waves, 
the respective far-from-equilibrium states are called op- 
tical turbulence [1-3]. For quantum physics, this is the 
simplest model of locally interacting bosons. The Hamil- 
tonian and the canonical NSE are as follows 

^ = /(|V^|^ + |^lV2)rfV, (1) 
= W'^n/5i}* = -A^) + IVI^^ . (2) 

To provide for a turbulent inverse cascade, we add 

medium-scale pumping and small-scale dissipation and 
solve the resulting equation numerically with a record 
spatial and temporal resolution. We address the central 
issue in the non-equilibrium physics, i.e. that of uni- 
versality: to what extent the properties of a far-from- 
equilibrium state depend on the excitation mechanism? 
Theoreticians prefer to pump a turbulent system by a 
random Gaussian force, which practically never can be 
met in real systems where waves usually appear as a 
result of some instability. Here we consider both ways 



to excite waves (additive random forcing and instability) 
and show that they produce similar turbulence only when 
the condensate is weak. We show that the instability- 
excited system undergoes a series of phase transitions 
to the states with different symmetries while the force- 
driven system docs not. 

Understanding the interaction of turbulence and a co- 
herent flow is an important problem in turbulence stud- 
ies in fluid mechanics and beyond, both from fundamen- 
tal and practical perspectives. In fluids, coherent flows 
are system-size vortices or zonal flows of diff'erent pro- 
files [4, 5], which are known to diminish the turbulence 
level, change its nature and make its statistics more non- 
Gaussian [4, 6]. Here we consider arguably the simplest 
case of a turbulent system with a condensate: the coher- 
ent part is expected to be a constant field and turbulence 
to consist of weakly interacting waves. We carry our sim- 
ulations well past time when all vortices, i.e. holes in the 
condensate, disappear (due to annihilations [7]) and ob- 
serve that the coherent part contains a small but impor- 
tant spatial structure and the turbulence is not effectively 
weak. As a result, this simple system demonstrates un- 
expectedly rich behavior with novel features never before 
observed in turbulence systems. 



II. THEORETICAL CONSIDERATIONS 

Let us denote = \tp\'^ and A^o = IV'P- what follows 
the overline means average over space, and the averages 
over time we shall denote by angular brackets. The total 
number of waves with nonzero wave numbers is N — Nq = 

The simplest condensate is a spatially-uniform field, 
A = ^/T^exp{—iNot), which is an exact solution of (2). 
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Small over-condensate fluctuations satisfy 

iijk = {k^ + 2No)i^k + A^r-k- (3) 
This equation gives the Bogolyubov dispersion relation, 

nl = 2iVofc' + , (4) 

for a pair of counter-propagating waves Vfe oc 
exp{—iNot — iflkt), V'* fe oc exp{iNot — iQkt)- 



A. Turbulence with the condensate 

As the condensate grows, the dispersion relation (4) 
approaches the acoustic one with the sound velocity 
c = \/2Nq. As known, acoustic waves running at the 
same direction interact strongly (producing shocks). On 
the other hand, the cfFcctivc matrix clement of the thrcc- 
phonon interaction behaves as |1423p — kxk2k^/NQ^^ 
(see e.g. [8]). One can estimate the mean value of the 
cubic term in the Hamiltonian using the weak-turbulence 
approximation [1]: 

ki,k2,k3 

=^ XI I^l23p'^l"2 <5(ki + k2 - k3)(5(f^i +D.2- ^3) 
\V\'^n^ck n^k 



k^ c A^Q^^ 



(5) 



Here wc; have introduced the mean spectral density of 
waves also called the normal correlation function rik = 
(IV'feP)- The factor k/c in (5) is the effective angle 
of interaction that can be estimated considering Q/. w 
c\k\ + k^/2c « ck+ck'^/2k+k^ /2c and comparing diffrac- 
tion and dispersion corrections: k±/k ~ k/c. We then 
estimate {'H2) — ckn and obtain the effective nonlinearity 
parameter for over-condensate phonons as the ratio: 



H2 



n 

Wo 



(6) 



This suggests that when the number of waves with zero 
momentum outgrows the number of waves with nonzero 
momentum, nonlinear interaction must be getting effec- 
tively weak despite the fact that the nonlinear term in 
the Hamiltonian is dominant. The estimate (6) explains 
the observation made in [9] that the statistics of the fluc- 
tuations is non-Gaussian at low condensate values but is 
getting closer to Gaussian, as the condensate grows. 

Formula (5) allows also to estimate the typical rate of 
a local nonlinear interaction (i.e. for waves with compa- 
rable k) as l/t„i{k) ~ \V\^nkk^{c/k^){k/c} ~ nuk/Nl'"^. 
The rate of nonlinear interaction l/i„; decreases as the 
condensate grows. Therefore, nonlocal interactions must 
play more and more important role. Here we define the 
coherent part as the field filtered at the frequency N and 



show that it is not a constant but has a rather elaborate 
spatial structure. Interaction of turbulence and the co- 
herent part (nonlocal in k-space) is responsible for the 
set of new phenomena described below. 

One may think that an effective weak nonlinearity ex- 
pressed by (6) allows one to describe the over-condensate 
turbulence by the weak turbulence theory i.e. by a closed 
kinetic equation written for the one-time (normal) cor- 
relation function rifc. We will see here that this is not 
necessarily the case. The crucial point is that validity of 
weak turbulence approach requires not only weak nonlin- 
earity but random phases as well. We show below (ana- 
lytically) that condensate imposes phase coherence, the 
most straightforward manifestation of which is the exis- 
tence of anomalous correlation functions and the phase 
coherence of the k, — fc-pairs with the condensate. We 
then describe the series of phase transitions (discovered 
by the direct numerical simulations) that show evolu- 
tion towards more and more ordered state as the con- 
densate grows: suc;h evolution cannot be described by a 
kinetic equation, which satisfies H-theorem and can de- 
scribe only entropy growth. 



B. Anomalous correlations 

Let us now show that an important element of tur- 
bulence against the background of a condensate must 
be anomalous correlations, i.e. phase coherence be- 
tween waves running in opposite directions (seems to 
have been overlooked so far). Indeed, the existence 
of the condensate must produce the anomalous cor- 
relation function {ipk'4'-k) oc ex.p{—2iNot) . To look 
into the effective dynamics of both the number of 
waves and the condensate, let us assume for simplicity 
that only two contra-propagating waves with amplitudes 
V'ifc = y/n±k exp{—iNot + i(f)±k) interact with the con- 
densate and expand the Hamiltonian (1) up to the terms 
quadratic in ip±k- That gives the equations: 



duk 



1 dN, 



dt 2 dt 
dB 

~dt 



= -B = 2NQ^nkn-k sin(2(/)o -(t>k- <i>-k) , 



4(2nfe + 2n-k - No ~ k^)C + SNoUkU-k 
-2N^{nk+n_k), 



dC 
Itt 



= (A^o + A;2 - 2nk - 2n_k)B 



where we have denoted A'^tplil^'^^, = C - iB/2. We 
see that the steady state corresponds to _B = and 
C = —Noy/nkU-k that is to 2(j)o — 4>k — 4>-k = tt (and 
Uk = ri-k = ^^/2)- Alternatively, one calculates the 
time derivative of the Hamiltonian, quadratic in A^ip±k: 
dU^/dt f» 4iVo/TO(A24^*^Vfe) = Ngy/HkU-k sin(2(/)o - 
(pk — (f>-k), which is the energy flux from the condensate 
to the wave pair with k, —k. The flux turns to zero and 
the steady state is established when the phase of the pair 
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is exactly opposite to the phase of A^, which plays the 
role of an effective pumping. 

C. Anisotropy of turbulence 

At the beginning, waves in the pumping shell are 
pumped isotropically. Eventually, the inverse cascade 
reaches lowest modes, which have the four-fold symmetry 
of the box. Let us now discuss what possible anisotropy 
may turbulence spectra acquire after the condensate ap- 
pears. Just looking at the Bogolyubov dispersion rela- 
tion (4) one may suggest the following scenario: spatial 
modulation of the condensate intensity A^o means respec- 
tive modulation of the sound velocity for shorter waves. 
Standing waves corresponding to the lowest modes of the 
box provide such modulation. Regions of low sound ve- 
locity would act as waveguides with short waves moving 
predominantly along the minima of A^o- That suggests 
that for sufBciently high Nq, when (4) is close to linear, 
the angular maxima at high k would correspond to min- 
ima at low k for a hypothetical anisotropic spectrum. 
In other words, acoustic waves interact effectively when 
they are coUinear, so the maximum of long waves at some 
angle effectively removes short waves from the pumping 
region at that angle. That consideration suggests a sim- 
ple picture of possible anisotropy: lowest modes have the 
symmetry of the square box and impose that symmetry 
on the small-scale turbulence (turned by 7r/4). The main 
result of this work is that this is not the case. 



III. NUMERICAL PROCEDURE 

Our numerical simulations are set up similarly to [9]. 
We use 4th order fully dealiased split-step method [10, 11] 
to solve Eq. (2) in the presence of forcing, 

iVt = -v^^-i-IV'I^V' + nV'- (7) 

The forcing consists of large scale pumping and small 

scale damping. The pumping-attenuation operator 7 is 
local in the k-space and is given by j{k) = jp{k) — jd{k) 
where the pumping is non-zero only within the shell [9], 

jp{k) = a ^ - A;f)(fc2 - k"^) at ki < k < kr 
and the damping is of Landau-type, 

jd{k)=/3k^h{k/kd), 
where function h{x) is defined as 

( (6x5)-i exp[5(l-x-2)] , 
h{x) = \ 

[ 1 - (5/6) exp [(1/2) (1 - x^)] , x>l. 

The pumping shell [fc; , /j^] is relatively narrow and occu- 
pies the middle of the computational domain. Note that 
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run B 
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8.00 X 10"=^ 


2.00 X 10"-^ 


run C 




1.60 X 10" 
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4.00 X 10"^ 


1.00 X 10"^ 


run D 




8.00 X IQ- 


-3 


2.00 X 10"^ 


5.00 X 10-* 



TABLE I: Note: (a) resolution is shown excluding dealiased 

modes; (b) the timestep is reduced to 5 x 10~' when max 
reaches 2 x 10^, and to 2 x 10"^ when max \%l>\'^ reaches 4 x 10 . 

we do not suppress long wave excitations (no damping is 
present at /c = 0) so, as expected, the condensate starts 
to form early on and the total number of waves begins 
to grow in time. As in [9], the initial condition is given 
by thermal equilibrium spectrum, rife = T/{k'^ + n), with 
random phases. 

The values for simulation parameters are shown in Ta- 
ble I. Note that the size of the box, L, increases with 
resolution, while forcing remains the same. As a result 
the number of waves per unit area, N, evolves in a similar 
way in simulations at different box sizes, see Figure la. 
Unless specified, the results are shown below for the run 
A. 

To collect better statistics, we occasionally stop the 
simulation and restart it with uniform friction which sta- 
bilizes the condensate growth and leads to a steady state. 
In this case, the forcing is replaced by 7 = 7 — 6, where b 
is constant determined empirically to stabilize A' at a 
given level. 

To study the structure of the coherent mode and fluc- 
tuations separately, we applied the time-spectral filtering 
used by Nazarenko and Onorato [7]. After the conden- 
sate appears, the frequency spectrum has a sharp peak 
at uo = N. We select a time window containing many 
27r/wo periods (typically 100) and store the temporal evo- 
lution of the ip{r,t) with 12-24 time-slices per period. 
This data is transformed to ui-k space, where the coher- 
ent mode — slice corresponding to coq — is separated. 
The rest of the frequency spectrum we treat as fluctua- 
tions. 

We use as a measure of nonlinearity, the dimensionless 
parameter e = where A is the pumping wave- 

length, which is not only the smallest scale in the inverse 
cascade but also the smallest scale in the system since 
dissipation dominates for the scales less than A. All our 
simulations are done with the same A « 27r/30, so e also 
has the meaning of rescaled number of waves, e « 0.01 A". 
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FIG. 1: Left: Number of waves as a function of time in L = 27r (dashed lines) and L = 47r (solid lines) simulations, also rescaled 
as e. Gray regions correspond to the transition between regimes; circled labels refer to the number of jets (petals). Right; 
Appearance of condensate in simulations with different box sizes in Run A. The condensate value A^o is shown by dots, next 
to the lines representing N. The initial linear growth of the condensate at strong forcing and large boxes is accompanied by 
creation and destruction of vortices. When vortices die out, No catches up with N. 



IV. RESULTS AND DISCUSSION 

A. Condensate growth 

The nonlinear interaction conserves the total number 
of waves so the latter grows only due to net pumping: 
dN{t)/dt = V{t) = J 'Ykfikit) (Pk. A simple expectation 
would be that N(t) must initially grow exponentially as 
a result of linear instability, then nonlinearity saturates 
this growth. After stabilizes in the pumping region, 
one expects a linear law TV = Ft. The first study of the 
wave number growth was done in [9] where it was ob- 
served that nonlinearity indeed slows it down (and for 
a while even reverses) after the initial exponential in- 
crease. The following assertions were then made in [9]: 
(a) the first nonlinear stage of growth follows the law 
N cx which means that the occupation numbers in 
the pumping shell decrease as cx (b) the square- 

root regime starts when the number of waves at the con- 
densate is of order of the total number of waves; (c) the 
square-root regime ends when e ~ 1 i.e. the correlation 
scale Nq ' is approaching the pumping scale; (d) after 
that, Uk stabilizes and N — Tt with time- independent 
r. As seen from the Figure 1 (left panel) (produced at 
the resolutions much higher and on time intervals much 
longer than those of [9]), the evolution, which we observe, 
does not follow this simple scenario. At earlier times, the 
growth on average can be crudely approximated as 
In some simulations, we observed a transitional linear 
regime, as in [9]. It can be seen on a linear scale, as, for 
instance, the section 500 < t < 1800 for L = Sir curve 
in Figure 1 (right panel). At later times, the growth 
rate of the condensate accelerates even faster than linear, 
with possible transition to the ultimate regime N oa t'^ 



regime hypothesized by Zakharov and Nazarenko [12]. 
The growth N{t) is determined by the number of waves 
in the pumping shell, and we shall see in the next sec- 
tion that this number evolves in a complicated way due 
to phase transitions. We found the law of this growth 
sensitive to different parameters and were not able to 
distinguish any definite scaling. That subject definitely 
requires extensive future studies. 

B. Phase transitions 

Initially, condensate is fed by an inverse cascade car- 
ried by an almost isotropic spectrum of fluctuations. At 
this stage, the nonlinear (interaction) term in the Hamil- 
tonian is less or comparable to the quadratic term (which 
describes dispersion for waves or kinetic energy for parti- 
cles). As condensate grows, the first symmetry breaking 
appears at e ~ 1 and the system becomes anisotropic 
— the amplitude and the phase develop patterns on the 
scale of the forcing and larger; the spectrum turns into 
an oval. The transformation happens gradually, in the 
interval 2 < e < 5, as the oval slowly becomes thinner at 
the waist. 

As the condensate amplitude grows further, we observe 
the sequence of phase transitions: at first the oval spec- 
trum turns into a dumbbell, then the symmetry changes 
from two-fold (two-petal) to three-fold and then to four- 
fold, as seen in Figure 2. Further transitions are possible. 
While we could not evolve our runs to higher levels of 
condensate, we observed other symmetries in the simula- 
tions with different initial conditions (discussed below). 
In Figure 1 (left panel) we have shown the regions of 
transitions as gray bands in the plot of e{t). The tran- 
sitions are accompanied by an increase in the variance 
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FIG. 2: Amplitude deviation from average in a 27r x 2n part of the computational domain (above) and spectrum (below) in 
L = Stt simulation. From left to right images correspond to t = 100, 1500, 2500, and 7500 (A^ = 219, 771, 1166, 4202, e = 2.40, 
8.45, 12.77, 46.08). 



transition 


two-fold to three-fold 


three-fold to four-fold 


run 


L 


t 


N 


e 


t 


N 


e 


A 


2-K 


2150 ±50 


932 ± 23 


10.22 ±0.25 


6250 ±50 


3910 ± 70 


42.88 ±0.77 




4-K 


2350 ±50 


1053 ± 24 


11.55 ±0.27 


7350 ± 50 


4310 ± 65 


47.30 ±0.71 




Stt 


2050 ± 50 


979 ± 26 


10.74 ±0.29 


7150 ±50 


3830 ± 57 


42.00 ± 0.63 


B 


2-K 
4-K 


13500 ± 500 
15250 ± 50 


1025 ± 45 
1010 ±5 


11.23 ±0.50 
11.08 ±0.05 


31500 ± 500 


2578 ± 79 


28.28 ±0.87 


N/A 


C 


2tt 


99500 ± 500 


1378 ± 7 


15.12 ±0.07 


N/A 



TABLE II: Time, the number of waves, and nonlinearity parameters at phase transitions. 



of the over-condensate excitations, n — N — Nq, as seen 
from the Figure 1 (right panel). 

The transition times and condensate amplitudes are 
listed in Table II for different pumping powers and box 
sizes. The transitions occur sharply yet the threshold 
values of e fluctuate from realization to realization, as we 
find in the simulations which only differ in the phases 
of initial thermal noise (see Figure 3.) This shows that 
we are not dealing with a linear instability with a well- 
defined threshold. Indeed, the condensate is stable with 
respect to the infinitesimal perturbations. It is likely 
that our transitions are of probabilistic nature similar 
to transitions in a pipe flow [13] or fiber laser [14] when 
with the change of a control parameter (condensate am- 
plitude in our case) one changes the probability that a 



finite-amplitude perturbation will lead the system away 
to a new state. 

It can be recognized in the second, third, and fourth 
panels of the second row in Figure 5 that the maxima at 
high k correspond to minima at low fc, in agreement with 
the waveguide argument presented in Section II C. How- 
ever, the symmetry of the spectra both at low and high 
k is definitely not that of a box. The minima and max- 
ima are directed at different angles and the spectrum of 
turbulence can have three-fold symmetry that is not that 
of a square box. Bottom row of Figure 4 shows that the 
mode coherently oscillating with the frequency N has an 
intricate spatial structure and multiscale correlations the 
anisotropy of which is directly related to the anisotropy of 
over-condensate fluctuations. The spectrum of the coher- 
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5000 




FIG. 3: The transition from 2 to 3 petals, seen as the kink 
in N{t), for different random phases of the initial condition. 
The insert plot is the zoom to the kink area. One can see 10% 
difference in A'^ at the moment of transition. 



ent mode has the doubled number of maxima and is gen- 
erally more symmetric than the whole spectrum, as seen 
from Figure 4. For nonzero momenta, over-condensate 
fluctuations dominate over the coherent part. 

In addition to phase transitions between different sym- 
metry states, the spectra undergo slow evolution between 
transitions. The intensity fluctuates, the whole structure 
may slowly rotate or swing. With time, two-petal spectra 
concentrate in a more narrow angle, as seen in Figure 6, 
and broaden again just before the transition. The three- 
petal spectra broaden with time, especially at high k. 
The third columns in Figures 4,5 show spectra just be- 
fore the transition, when they become close to six-petal; 
the seeds of the new (four-fold) symmetry visible at low k 
in the coherent part, as seen in Figure 7. The population 
of the pumping shell increases between the transitions 
and falls during the transitions, which explains a compli- 
cated character of the evolution of the total number of 
waves. 

At any moment, except possibly very close to the tran- 
sition events, and for all symmetries, the angular width 
of the spectra decreases towards larger wavenumbers, as 
seen in Figures 4, 5. This may be analogous to spectrum 
anisotropization observed for acoustic turbulence [1, 15]. 
It is worth stressing the difference in two cases: acoustic 
turbulence studied there corresponds to a direct energy 
cascade (realized by pumping-generated long waves turn- 
ing into shocks), while here we have an inverse cascade. 
One can therefore say that as an inverse cascade pro- 
ceeds towards larger wavelength, the spectra are getting 
wider as weak turbulence theory would predict [1]; yet 
presently we see no way it can predict the spontaneous 
appearance of anisotropy and changes in symmetry, par- 
ticularly, signiflcant anisotropy at the pumping scale. 

One may suggest a possible physical mechanism of 



phase transitions as follows: phonons effectively inter- 
act only within the angle k/^/No, which decreases as Nq 
grows, so the turbulence tends to be broken into jets. 
The number of jets j increases proportionally to the in- 
verse of the interaction angle i.e. as iVp/^ cx Indeed, 
we observe that the transitions happen approximately at 
Cj ~ j^. Even if this is indeed the basic mechanism of the 
transitions, we still lack any understanding of whether 
there is a unifying principle that can predict which tur- 
bulent state is realized the way variational principle does 
it for thermal equilibrium. One direction worth explor- 
ing is whether one can develop an approach similar to 
the weak crystallization theory [16] despite the fact the 
we are dealing with developed turbulence and power-law 
spectra carrying flux in k-space, as seen in Figs 5, 11 
below. 

Finally, we made separate runs taking initial conditions 
for -0 as a thermal noise plus large constant (pre-existing 
condensate). In these cases, the same states (with 2,3,4 
petals) appear i.e. they are true attractors, independent 
of the history. We also observed 6-petal spectra with 
triangular spatial pattern in the amplitude, and 8-petal 
spectra with the spatial distribution comprised of patches 
of square patterns oriented at 45 degree angle. The 6- 
petal spectra was observed at higher e than 4-petal spec- 
tra, although pre-existing condensate makes all thresh- 
olds lower and possibly less sensitive to the phases of the 
initial thermal noise. We speculate that this is because 
in the runs with a pre-existing condensate, the spectrum 
first develops a uniformly populated pumping shell which 
is then free to break apart into the appropriate number 
of jets, contrary to evolving systems where the previous 
configuration has to re-populate the pumping shell. The 
ability of the system to quickly find its equilibrium state 
opens horizons for a much more extensive exploration 
of parameter space, which hopefully will be a subject of 
future studies. 



C. Orientational order 

For two-fold and three-fold spectra, turbulence statis- 
tics is translational invariant. The most remarkable find- 
ing is seen at the last image of the top row in Figure 2: 
on top of a long-range orientational order, a short-range 
positional order appears. This state is between the solid 
and the isotropic liquid very much like hexatic phase in 
2d melting [17]. It requires future detailed studies to es- 
tablish whether this is indeed a turbulent analog of the 
Berezinski-Kosterlitz-Thouless transition; to avoid mis- 
understanding, let us stress that there are no vortices 
(holes) in the condensate at this stage, the phase fluc- 
tuates weakly; it is the small condensate perturbations 
which are getting ordered. These spatial correlations are 
also seen in Figure 8 which shows the correlation func- 

2 

tion of the amplitude, %li{x)ip{x + r) — ilj , where overline 
means averaging over the position x. 

The type of crystallization that we observe is very 
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FIG. 4: Total spectrum (above) and filtered spectrum of the coherent mode (below) in L = Stt simulation. From left to right 
images correspond to t = 1500, 2100, 6500 and 7500 {N = 771, 1005, 3102, 4202). The averaging was done over 100 time units 
in stabilized runs. Colormap is the same as in Figure 2. 




Ji/2 71 3k/2 2it 7c/2 71 37t/2 2jt it/2 k 3%i2 2% Jt/2 % 3n/2 2% 

angle angle angle angle 

FIG. 5; Radial dependency at two directions (first row) and angular dependency at |fc| = 8, 24, 48, 64 (second row, from top 
to bottom) for the spectra shown in Figure 4. Dashed lines in the radial spectra correspond to the coherent mode. 



much different from the vortex crystallization with a D. Anomalous correlations and collective 

wavelength externally imposed by a cut-off in 2d incom- oscillations 
pressible turbulence [18]. 

As we have predicted above (Section II B) anoma- 
lous correlations must exist between any two contra- 
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FIG. 6: Form of the spectrum at t = 10000, 40000, 80000, 
from left to right, in Run C at L = 27r. The spectrum angular 
width decreases as grows. 




FIG. 7: Zoom into the central part of the total spectra (above) 
and the spectra for the coherent mode (below) for L = Att and 
t = 3000, 6000; both correspond to three-fold symmetry. One 
can see the four-fold symmetry appearing at low k in the lower 
right panel. 



propagating waves. In our direct numerical simulations 
we have confirmed that this is true at least for the two 
lowest modes of the box, their total phase is indeed 
shifted by tt relative to the condensate. These modes 
give the main contribution of the collective oscillations of 
the system which are periodic conversions of waves from 
condensate to over-condensate with the total number of 
waves unchanged. 

We checked that in the numerical data the steady- 
state anomalous correlation function indeed has the 
phase shifted by tt relative to the condensate value 
squared. For example, for L — Att and N = 570, the 
modes that have the largest amplitudes, with {kx,ky) — 
(1, 0), (1, 1), (1, — 1), have their phase differences {(j)k + 
(t>-k - 20o)/7r = 1.01, 1.02,0.99 respectively 

It was noticed in [19] that when broad turbulent spec- 
tra coexist with a sharp spectral peak these two groups 
of waves generally have different relaxation times which 
makes possible oscillations about a steady state. For 
three-wave interaction, an integral model (of a predator- 
prey type) that describes the evolution of the total num- 
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FIG. 8: Correlation function of the amplitude of the over- 
condensate fluctuations at f = 7500 for L — 2n. 

bers of waves in the two groups has the form: dNo/dt = 
—bNo + Nqu, dn/dt = 771 — Nqti, which gives the oscilla- 
tions with the frequency \/jb — \J uNq [19]. 

The oscillations for NSE with a condensate were re- 
ported in [8] where a similar model was suggested: 
dNo/dt = -bNo + Nqh'^, dn/dt = 771 - N^n. This model 
gives the steady-state values = .^7, fi = ^/b and 

the same frequency of oscillations Vt^ — \/ uNq. The 
most evident defect of the latter model is that it does 
not conserve the total number of waves, moreover, the 
steady-state values are wrong too (in reality, Nq oc 1/6 
while n is practically independent of b in wide intervals) . 
We notice here that to describe the collective oscillations 
of the turbulence-condensate system, one needs to ac- 
count properly for the anomalous correlation function, 
as was done while deriving system (7). In particular, our 
suggested system does describe the oscillations around a 
steady state. Consider, A^o = ^0 + Uk — fik + m, 
assuming that the steady-state values satisfy Nq ^ fik- 
We then get 

dm _ ^_ dM 
~dt ^ 2dt ' 

^^iNo + k')B 

The system above describes collective oscillations 
with twice the Bogolyubov frequency: M^m, B,C oc 
exp(-im), n ^2nk ^ 2V2iVofc2+l!4. Indeed, numer- 
ical simulations show that the oscillations have frequen- 
cies that are equal to twice that of the lowest modes of 
the Bogolyubov spectrum. For the box with L — 2tt 
we find that the modulus of the anomalous correlation 
functions a{kx,ky) = cr(l, 0) and cr(l, 1) oscillate respec- 
tively with the frequencies 2^/2Nq and A^/JT^. These 
frequencies (which are much smaller than the frequency 
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FIG. 9: Left: oscillations in over-condensate in Run A, L = 2ti box, zoomed around times t' = 800 (A'' = 486) and t' = 3500 
{N = 1518). Middle: frequency spectrum of the time-dependence of the number of over-condensate waves for the same run in 
time intervals [800, 810] and [3500, 3510]. Right: first two frequencies as function of the number of waves in L = 2ii and L = 47r 
boxes. Dashed lines have the slopes 4, 2\/2, 2, \/2- 



6090 



6080 



6070 



6040 Pi: 



6030 



coherent mode, all k 




coherent mode, k = 



0.05 



0.1 

t-t' 



0.15 



0.2 



FIG. 10: Oscillations of number of waves in the condensate. 
No, and in the coherent mode for L — 2n. 



of the condensate N) essentially do not depend on the 
level of over-condensate fluctuations, they are clearly 
seen in the oscillations of the condensate amplitude and 
of the normal correlation functions shown in Figure 9. 
As seen from Figure 10, the total number of waves 
is unchanged. The number of over condensate waves 
(about 40 — 6080 — 6040) consists of a coherent part 
(about 10 = 6050 — 6040), the rest (about 30) are waves 
with other frequencies. Note that the number of over- 
condensate waves exceeds substantially the amplitude of 
the oscillations (about 5). The oscillations can be thus 
treated in a linearized approximation. 



E. Moments, spectra and flux lines 

Table III shows the integral statistical characteristics of 
over-condensate fluctuations for different times and levels 



t 


N 


No 


\5i>\yNo 


\5^^\y[\& 


1500 


771 


751 


0.027 


2.95 


2500 


1166 


1133 


0.029 


2.67 


6500 


3102 


2871 


0.080 


3.17 


7500 


4200 


3986 


0.054 


2.40 



TABLE III: Integral statistics of fluctuations corresponding 
to spectra in Figures 4,5. 



of the condensate. We define Sip = ip — iJj and compute 
the spatial averages of its moments. The fourth column 
gives the second moment, which is the fraction of waves 
with nonzero wavenumbers: \6ip\'^/No = (TV — Nq)/No. 
The last column shows the flatness, which is close to the 
Gaussian value 3. Far from the transitions, the flatness 
is actually somewhat lower than 3 and decreases as the 
condensate grows, i.e. the condensate effectively sup- 
presses strong over-condensate fluctuations contrary to 
what condensate vortices do to 2D incompressible tur- 
bulence [4]. Only for t = 6500, which is close to the 
three-to-four transition, we find flatness slightly exceed- 
ing 3, i.e. somewhat stronger fluctuations. Note in pass- 
ing that the coherent mode at all scales has substantially 
smaller flatness, which means that strong fluctuations of 
the mode are even more rare than for a Gaussian random 
quantity. 

Let us stress that despite anomalous correlations, 
symmetry-breaking phase transitions and collective os- 
cillations (all properties normally associated with narrow 
spectra in k-space) , our spectra are wide along and across 
the jets, decaying approximately by power laws, as seen 
in Figure 5. Note that is the spectrum of weak tur- 
bulence inverse cascade which coincides in this case (up 
to a logarithmic factor) with the thermal equilibrium, 
while is the spectrum due to the shock waves. 

Let us now describe the lines of the flux that these spec- 
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FIG. 11: Source distribution and flux flow lines a,t t — 1500, 
t = 2500, and t = 7500 for L = Stt. 




FIG. 12: Spectra for an instability (left) and the random force 
(right) at about the same levels of condensate for L = 2ti. One 
sees that force-driven turbulence remains isotropic except for 
a few box-size modes, while the instability-driven turbulence 
develops a three-petal spectrum in the whole of k-space. 

tra carry towards the condensate. Conservation of waves 
allows one to write NSE as a continuity equation and de- 
fine the flux, Jfc, via the relation divJ^ = (|*AfcP)7fc = 
Qk- In this Poisson equation, Qk plays the role of a 
charge density and the flux plays the role of potential. 
We compute Qk and solve the above equation. Figure 11 
shows the flux lines while positive/negative source Qk is 
shown by red/blue regions. 

F. Non-universality of the inverse cascade 
turbulence 

It is important that the phase transitions happen only 
when the waves are excited by an instability. If we ex- 
cite waves by an additive random force concentrated in 
the same ring in k-space with about the same amplitude, 
turbulence remains close to isotropic at all wavenumbers 



exceeding unity and at all condensate levels, as seen in 
Fig 12. To the best of our knowledge, this is the first 
demonstration of significant non-universality of turbu- 
lence with respect to excitation mechanism. It is par- 
ticularly striking that we find it for an inverse cascade 
which is generally particularly robust [20]. Apparently, 
the condensate-turbulence interaction can make turbu- 
lence non-universal. 

This set of simulations was performed for L — 2t: and 
the following procedure was used. First a run with an in- 
stability excitation was done with the initial condition as 
a large constant (corresponding to Nq = 1000) plus small 
thermal noise. The run was stopped after it produced a 
3-petal spectrum shown in the left panel of Fig. 12; dur- 
ing that time the total number of waves grew only by 
3.6% from its initial value. Second, a run with an addi- 
tive forcing was performed. In this run the multiplicative 
term 7p(A:)'0fc was replaced by an additive forcing /(k, i) 
(the multiplicative damping term, jd{k), was preserved). 
This forcing term was non-zero only inside the pumping 
shell ki < k < kr, it had random phases uniformly dis- 
tributed in the interval [0, 27r] and uncorrclated both in 
k-space and time. Its amplitude was chosen to be con- 
stant |/(k, t)) \ — f the value of which was determined by 
a condition that both additive and multiplicative forcing 
have the same average amplitude in the k-space: 

•^"^TP^^ / lp{k)\^Pk\d^k . 

-KlKj, — ) Jki<k<k^ 

To conclude, discerning guiding principles that govern 
the physics far from equilibrium remain elusive. Is there 
any quantity (similar to entropy) that the system tries 
to optimize by undergoing the series of phase transitions 
described here? 
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